On the Matrix Representation of Quantum Operations 
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This paper considers two frequently used matrix representations — what we call the \~ an d S- 
matrices — of a quantum operation and their applications. The matrices defined with respect to an 
arbitrary operator basis, that is, the orthonormal basis for the space of linear operators on the state 
space are considered for a general operation acting on a single or two d-level quantum system (qudit). 
We show that the two matrices are given by the expansion coefficients of the Liouville superoperator 
as well as the associated bijective, positive operator on the doubled-space defined with respect to two 
types of induced operator basis having different tensor product structures, i.e., Kronecker products of 
the relevant operator basis and dyadic products of the associated bipartite state basis. The explicit 
conversion formulas between the two matrices are established as a computable matrix multiplication. 
Extention to more qudits case is trivial. Several applications of these matrices and the conversion 
formulas in quantum information science and technology are presented. 

PACS numbers: 03.67.-a, 03.65.Wj 



I. INTRODUCTION 

The formalism of quantum operations offers us a pow- 
erful tool for describing the dynamics of quantum sys- 
tems occurring in quantum computation, including uni- 
tary evolution as well as non- unitary evolution ElQ- It 
can deal with several central issues: measurements in the 
middle of the computation, decoherence and noise, us- 
ing probabilistic subroutines, etc. 0- It describes the 
most general transformation allowed by quantum me- 
chanics for an initially isolated quantum system 0, 
Experimental characterization and analysis of quantum 
operations is an essential component of on-going ef- 
forts to develop devices capable of reliable quantum 
computing and quantum communications, and is a re- 
search subject of considerable recent interest. There have 
been extensive efforts on the statistical estimation of 
quantum operations occurring in natural or engineered 
quantum processes from experimental data known as 
"quantum channel identification" or "quantum process 
tomography" @ S i [H [H [12 . 

There are several ways of introducing the notion of 
a quantum operation, one of which is to consider it as 
superoperators acting on the space of linear operators 
0, ll'il ll l| . Any physical quantum operation has to be 
described by a superoperator that is completely positive 
(CP); i.e., it should map the set of density operators 
acting on the trivially extended Hilbert space to itself 
0, • It is known that any CP map can be decomposed 
into the so-called Kraus form by a set of Kraus opera- 
tors 0. However, this description is unique only up to 
unitary equivalence 0, just like the decomposition of a 
given density operator into convex sum of distinct, but 
not necessarily orthogonal, projectors is unique only up 
to unitary equivalence 0, Il5| . Alternatively, the super- 
operators can be represented in matrix form by provid- 
ing an operator basis i.e., the orthonormal basis for 
the space of linear operators on the state space, just as 



the operators on the Hilbert space can be represented in 
matrix form by providing a state basis for the Hilbert 
space. For example, the density operator can be repre- 
sented by a density matrix defined with respect to the 
chosen state basis. The density matrix provides a unique 
description of the quantum state once the state basis is 
fixed, although we still have a freedom in choosing the 
state basis. Similarly, a quantum operation can also be 
uniquely described using the matrix once the operator 
basis has been fixed. 

One can construct many different types of matrix rep- 
resentation. This paper considers two different matrix 
representations for superoperators frequently found in 
the literature. The first one is what we call the ^-matrix, 
which is also called the process or dynamical matrix by 
several authors [1 Efl El El El III Let us con- 
sider a d-dimensional Hilbert space (-ff-space) 7i<j, and 
the space of linear operators acting on Tid with a scalar 
product (A,B) = Tr A*B, that is, the Hilbert-Schmidt 

space (HS-spa.ce) TLSd- If we choose the fixed basis set 
.2 i 

{E a } a=0 in HSd, the linear operation 6 can be repre- 
sented by the binary form of the superoperator 



^(©) = Xa ^ a ®l 



(i) 



a,0=O 



acting on TLSd, which maps a linear operator in TLSd 
into another one. In Eq. ijT)!. the substitution symbol 
© should be replaced by a transformed operator, and 
double-hat * is used to distinguish the superoperator from 
an ordinary operator acting on Tid- The coefficients Xaf3 

d 2 1 

form a d 2 x d 2 positive matrix \ = [Xa/3] a p =0 , if 6 is a 
physical quantum operation plj . 

Alternatively, another matrix representation of a quan- 
tum operation is given in terms of the Liouville formal- 
ism |23, |23| • In this formalism, the linear operators 
in TLSd are identified with the supervectors in a Liou- 
ville space (L-space) Ld? . Introducing a double bra-ket 



2 



notation for the elements of Cd 2 [23j , we associate every 
operator A with an L-ket \A)) and its Hermitian con- 
jugate operator A 1 with an L-bra ((A\. The space L d i 
is furnished with an inner product ((A\B)) = Tr A^B, 
and constitutes a d 2 -dimensional Hilbert space. Then by 
chosing an arbitrary fixed set of operator basis {E a } a=0 
in TCSd, any linear operation 6 can be written as the 
superoperator 

d 2 -i 

J= J2 S a p\E a )){{Ep\ (2) 

a, 0=0 

acting on C,p , which maps a L-space supervector into 
another one. The coefficients S a p form a d 2 xd 2 complex 

matrix S = [»5 a| a]^ p= - They represent the amplitudes 

of the operator components E a contained in the state 
after applying the quantum operation on the operator 
component Ep. We call it the iS-matrix by analogy with 
the S-matrix appearing in time-independent scattering 
theory 0, 113 • The 5-matrix has been actually used 
to describe quantum operations by several researchers 
on quantum information science H, H |2(| . The %- and 
5-matrices offer us the most general description of the 
dynamics of initially isolated quantum systems allowed in 
quantum mechanics, just as the density matrix offers us 
the most general description of the quantum mechanical 
state. 

The choice of matrix representation type is a matter of 
convenience, depending on the application. We will dis- 
cuss later how the x~ arL d 5-matrices are useful for the 
analysis and design of quantum operations. Although 
these matrices indeed have their own useful applications, 
their mutual relation is non-trivial and has not been clar- 
ified. The main purpose of this paper is to clarify the 
underlying relation between two different matrix repre- 
sentations of quantum operations, and to provide the way 
for building bridges across the different classes of appli- 
cations. We here consider a quantum operation acting 
on the state of a single d-level quantum system (abbre- 
viated as single-qudit operation) or a two d-level quan- 
tum systems (two-qudit operation). We start in Sec. |H] 
by recalling the notion of operator basis and its proper- 
ties, which is helpful for the subsequent discussions. We 
note the equivalence between the supervectors in the L- 
space Cfp and the vectors in the doubled Hilbert space 
r H® 2 =Hd H-d- This equivalence implies that for any 

,2 -. 

operator-basis set {E a }^~ for TtSd, there is an isomor- 
phic state-basis set {lEa}}}^'^ iorHf 2 . We review sev- 
eral properties of the operator basis for later discussion. 
In Sec. IIIII we first consider the single-qudit operations. 
We show that x~ an d 5-matrices are given by the ex- 
pansion coefficients of the L-space superoperator S and 

the associated operator \ = «5 ®i(dpi) acting on Tif 2 
defined with respect to two types of induced operator 
basis on 7~t® 2 having different tensor product structures. 
Here, pi is the density operator of the isotropic state 



on Hf , and 1(0) = is the identity superoperator 
acting on the HS-sp&ce of the second system. This re- 
sult implies that there is a bijection between S and x, 
from which we can deduce the conversion formula be- 
tween x~ an d iS-matrices as a computable matrix alge- 
bra. Although Nielsen and Chuang have considered such 
a formula [2 IT^ |. their method requires finding matrix 
inverses to convert from the 5-matrix to the ^-matrix. 
Here, we show a conversion formula without matrix in- 
version. We then extend the formula to two-qudit op- 
erations. We also briefly review the requirement for the 
X-matrix to represent physical operations. In Sec. II VI 
we illustrate the applications of the present formulation. 
First, we discuss how x~ an d 5-matrices can be obtained 
experimentally. We describe a typical procedure to ob- 
tain the x~ an d iS-matrices defined with respect to an 
arbitrary operator basis set. Next, we discuss how these 
matrices and the present conversion formulas are useful 
for the analysis and design of the quantum operations, 
quantum circuits, as well as quantum algorithms. In Sec. 
Ivl we summarize the results. 



II. OPERATOR BASIS 

It has been noted by many researchers that the su- 
pervectors defined in the L-space £^2 can be identi- 
fied with the vectors in the doubled Hilbert space Ttf 2 
E3, IH |M HJ. Let us start by reviewing this fact, 
briefly. Consider an arbitrary chosen set of an orthonor- 
mal basis {l^)}^ 1 for Hd (denoted as standard state ba- 
sis). Any linear operator A in TLSd can be expanded 

as A = Y^jLo-Aijlfytf] where the dyadic products 
form a basis for TLSd- The L-space supervectors corre- 
sponding to the dyadic operator is denoted by the 
double ket \ij)) , with which the L-space supervector as- 
sociated with A is written as \A)) = X}j7=o ^ylij))- The 
scalar product of two L-space supervectors \A)) and \B)) 
is defined as 

- Trits, ( 3 ) 

which introduces a metric of an L-space. The vectors 
form a basis for a ci 2 -dimensional Hilbert space. Thus, we 
can safely identify with the product state |i) \ j). 
Then, the L-space vector associated with A can be identi- 
fied with vector \A)) = (A I)\ I)) in the doubled space 
Hf IH, where d" 1 / 2 !/)) = <r 1/2 Ei=o |*> |») is 
the isotropic state in Tt® 2 [H H3 , a nd 7 is the identity 
operator in TLSd- It may be helpful to recall the mathe- 
matical representations of A in the space C dxd of d x d 
complex matrices and the vector \A)) in C d . Consider a 
representation where \i) is a column vector in C d with a 
unit element in the jth row and zeros elsewhere. Then, A 
is identified with the dx d complex matrix A = [Ay]^ 

in C dxd , and \A)) is obtained by placing the entries of a 
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d x d matrix into a column vector of size d 2 row-by-row, 
i.e., 

\A)) = [A u ,--- ,A ld ,A 2U --- ,A 2d ,--- ,A dl ,--- ,A dd ] T . 

(4) 

Therefore, |^4)) contains the same elements as A but in 
different positions. This and Eq. @ indicate that A and 
\A)) are isometrically isomorphic. Accordingly, we may 
identify the Z^space with the doubled Hilbert space, i.e., 
C d 2= H® 2 . Hereafter, we use a common symbol T~t® 2 to 
denote both these spaces without loss of clarity. 

It will be useful for the later discussion to note the 
following relations hold: 

A®B\C)) = \ACB T )), (5) 



Tr 2 [|i)> 12l2 <<S|] = (i£t)M (6) 



Tr 1 [\A)) 12l2 ((B\] = (A T B*) (2 \ (7) 

where the indices refer to the factors in Tif 2 in which the 
corresponding operators have a nontrivial action, and the 
transposition and co njug ation are referred to the chosen 
standard state basis |16| . 

Now, let us consider the operator basis, that is, the 
complete basis for HS d . Consider an arbitrary set of ca- 
veators in Ti.f 2 . From the above isomorphism, this set 
can be written as {\E a }}}^S~J , where {E a }^~g is the 

associated set in 7iS d . The set {\E a ))}^S^ is the state 
basis for Tif 2 iff it is orthonormal, i.e., 

((E a \E p )) = 5 aP , (8) 
and it is complete, i.e., 
d 2 -i 

J2\E a ))((E a \=I®I. (9) 

The previous discussion shows that the vector \A)) in 
Tif 2 is identified with T^space supervectors associated 
with the operator A in 7iS d , This implies the following 
proposition. 

,2 i 

Proposition 1 A set of the operators {E a } a=0 is a ba- 

A ,2 -I 

sis setforHS d iff a set of states {\E a ))} a ~ is the basis 
set for Hf 2 . 

We argue below that this is true. First, we note the 
following lemmmas. 

Lemma 1 A set of the operators {i? a } Q=0 in 7iS d is 

orthonormal iff a set of states {{Ea))}^^ in 7i. d 2 is 
orthonormal. 



This is a trivial consequence of Eq. (J3J) . 

Lemma 2 A set of the operators {E a } a=0 in 7iS d is 

complete iff a set of states {\E a ))}^=Q in Ti. d 2 is com- 
plete. 

To prove this, the following Theorem is helpful [Tt| : 
Theorem 1 (DAriano, Presti, and Sacchi (2000)) A 

,2 -I 

set of the operators {£" Q } Q=0 in 7iS d is complete iff it 
satisfies one of the following equivalent statements: 

1. For any linear operator A on H. d , we have 

d 2 ~i 

i=£(Tr£ti)£ a . (10) 

2. Let £ de poi{' • • ) be the superoperator on the space of 
linear operators in TLS d describing completely de- 
polarizing operation. For any linear operator A on 
7i d , we have 

1 1 

£de P oi(A) = ^(Tri)/ = - E a AEl (11) 

3. For chosen any state basis {\i)}iZi forH d , we have 

d 2 -i 

(n\El\m)(l\E a \k) = S nk S ml . (12) 

Now, let us prove Lemma|5] If {l-E'a))}^^ 1 i s complete, 
Eq. 10 must be satisfied. Then, for any A in HS d , we 
have 

d 2 -i 

(A®i)\i)) = \A)) = J2 \E a ))((E a \A)) 
d 2 -l 

= Y(TrElA)E a ®i\i)). 

Since this holds for any A, Eq. (|10fl must be satisfied. 
Hence, {E a }^~Q is complete. Conversely, if {E a }^~Q is 
complete, then for any \A)) in Ttf 2 , we have 

d 2 -i 

\A)) = (A®I)\I)) = £(Tr£ti)£: a ®/|7>) 

a=0 
d 2 -l 

= \E a ))((E a \A)). 

a=0 

Since this holds for any \A)) , Eq. © must be satisfied. 

,2 -J 

Hence, {\E a ))} a=0 is complete. 
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From these two lemmas, we obtain Proposition^] This 
indicates that the state basis \E a )) for Ttf 2 has bijective 

correspondence to the operator basis E a for HSd. The 
following corollary is a consequence of Theorem ^ which 
will be useful for the later discussions. 

,2 -I 

Corollary 1 Let {£ a } Q=0 be an arbitrary chosen set of 
an operator basis for TLSd ■ The isotropic state in Ttf 2 is 
written as 

p I = ~\i))((i\ = ^E a ®E*, (13) 

ct=0 

and the swap operator Von T~if 2 is written as 

V=Y,E a ®El (14) 

Equation i|13|) can be proven by explicit evaluation of 
the matrix elements and using Eq. (|12|l . Equation i|14|) 
is obtained by performing the partial transpose on both 
sides of Eq. (|13fl with respect to the second system. 

Examples of the operator basis 

We show three illustrative examples of the operator 
basis for Tid frequently found in the literature. The 
first example is a set of transition operators Ttuj) ■— 
\i)(j\ with i,j = 0, - • - , d — 1, and (i,j) := di + j 

[34|. The associated states form a basis set {\^ a )}}a=o 
whose elements are the tensor product of the standard 
state basis, i.e., = \i) <8> \j). The next ex- 

ample is a set of unitary irreducible representations of 
the group SU(d) or the discrete displacement operators 
on the phase-space torus, whose elements are J7( TOi „) = 
w mn/ 2 u mk t (mn , k) /Vd with m,n = 0, • • • , d - 1, 
where u> = \ x / d = e l2 ^^ d and © denotes addition mod- 
ulo d |3 HE Hll • The associated states form a basis set 
{\U a ))}a=o of G? 2 -orthogonal maximally entangled states 
in Ttf 2 . The last example is a set of d 2 ~ 1 traceless Her- 
mitian generators of the group SU (d) supplemented with 

,2 -. 

the normalized identity operator given by {A Q } Q=0 = 

{I/Vd, Uo,l,MO,2,- ' ' ,Md-2,d-l,«0,l, #0,2, "' ' ^d-2,d-l, 

Wi,W2,--- ,Wd-i} where d(d — 1) off-diagonal genera- 
tors are given by u it j = (tt^j) + K(j^))/\/2, v it j = 

— Tt(j,i))/V2 with 0<i<j<d — 1, 
and d — 1 diagonal generators are given by w k — 
(- y^ -n 7T(i,i) + k% {k>k) )/ \/k{k + 1) with 1 < k < d - 1 
|34| . The choice of basis is of course a matter of conve- 
nience, depending on the application. 

Since the associated sets {\Tr a ))}t=o > {l^a))}^^ 1 ' an d 
{\^a))}a=o ar e state bases for Tif 2 , they should be uni- 
tarily related. This implies that the operator bases 7r Q , 
U a , and Aq, should also be unitarily related. In general, 
two sets of state basis {I-Eq))}^ 1 and {\F a ))} d ~^ in 



Tlf are unitarily related, i.e., 

d 2 -l d 2 -l 

\F P )) = Y, \E a )){{E a \F (j )) = Y \E a ))U aP (15) 

iff the operator bases E a and F a in HSd are unitarily 
related, i.e., 

Fp=J2 E a U a0 . (16) 

In Eqs. JTSJ and JTSJl, W Q/3 = ({E a \Fp)) = TrElFp is a 
a/3-entry of the d 2 x d 2 unitary matrix U. If we consider 
a unitary superoperator acting on the vectors in Hf 2 , 

d 2 ~\ 

U= E U a p\E a ))((E \, (17) 

q,/3=0 

Equations. I|15|) and (|16|l are the unitary transformation 
of the operators. Note that Eq. Q16[l does not imply uni- 
tary equivalence of E a and F a , i.e., Fp = WE a W' for 
some unitary operator W in HSd, although the unitary 
equivalence of E a and F a implies Eq. (|16|) . In general, 
F/3 7^ WE a W^ for any unitary operator VF in HSd, even 
if Eq. I|16|) holds. For example, 7T(i.j) is the operator ba- 
sis of all the elements of which have rank one, whereas 
U a and A a are those operator bases of all the elements 
of which have rank exceeding one. Therefore, TT(ij) is in- 
duced by a standard state basis {l*)}^ 1 for Tid, whereas 
U a and A a are not. This clearly indicates that either the 
set {U a }a=o or {AqIq^ 1 i s unitarily related to the set 

d 2 — 1 

{■^a} Q= o ' but is not unitarily equivalent to this set. 

III. MATRIX REPRESENTATION OF 
QUANTUM OPERATIONS 

Let us turn our attention to a single-qudit operation. 
As shown in Sec. [IJ this quantum operation can be rep- 
resented by either an ffi>-space superoperator or an L- 
space superoperator, in which the \~ an d 5-matrices are 
introduced with respect to the arbitrary, but associated 
sets of basis {E a } d ^ for HS d and {\E a ))} d ^ for H® 2 , 
respectively. In this section, the underlying relationship 
between these two matrix representations is discussed. 

Let us first consider the L-space superoperator. In the 
Liouvillc formalism, the operators in HSd are identified 
with the vectors in Ttf 2 . Any operation & is identified 

with the one-sided operator S acting on Hf 2 , which can 
be expanded using the state basis \E a )) for Tif 2 as shown 
in Eq. @. The elements of a d 2 x d 2 matrix S are for- 
mally written as S a p = ((E a \S\Ep)) . Alternatively, the 
same operation 6 is written as a two-sided superoperator 
acting on the operator in HS d , which can be expanded 
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using the operator basis E a for TLSd as shown in Eq. . 

Since S\Ep)) — \S(Ep))), we find the matrix element S a p 
can also be written as 

d 2 -i 

7,<5=0 

where we used Eq. (JSJ. Substituting the right-hand side 

of Eq. JTSJ for S a)3 in Eq. ©, we find that S can be 
written in terms of either the matrix 5 or % as 

d 2 -l d 2 -l 

£= S ^\Ea))({E \= XapE a ®E* . (19) 

a, 13=0 a, (3=0 

In Eq. (|19|l . we find two types of induced operator ba- 
sis on Hf 2 having different tensor product structures, 
that is, Kronecker products E a <g)Ep and dyadic products 
\E a ))((Ep 1 of the state basis associated with the operator 

basis set {E a } a=0 . Note that both types of basis set do 
not cover all the possible basis sets on TL® 2 . Obviously, 
the former type of basis set covers only those sets that 
are factorable with respect to the original and extended 
system spaces. For example, the set of d 4 -dyadic prod- 
ucts of the <i 2 -maximally-entangled states in Tif 2 is not 
a factorable basis set, and can not be covered by the for- 
mer type. Similarly, only an operator basis on Tif 2 with 
all elements of which have rank one can be reduced to 
the latter type of the basis set. Therefore, each type of 
basis set can describe its own particular subset of all the 
possible basis sets on TL® 2 . 

Let us next consider another operator x on Ttf 2 , which 
we call the Choi operator j^. We will show that this 
operator has bijective correspondence to the L-space su- 

peroperator S. It is known that the isomorphism be- 
tween the operator in TLSd and the bipartite vector in 
TL® can be straightforwardly extended to the isomor- 
phism between the superoperator acting on TLSd and the 
operator acting on Ttf 2 . Jamiolkowski first showed that 

the map between the HS-sp&ce superoperator S (©) and 
the operator \ = <5 © T{dpi) acting on TLf 2 is an iso- 
morphism, where S (©) and X (©) act on the i/5-space 
of the first and second systems, respectively j38(. If we 
note Eq. i|13fl an d the following equivalent relation that 
follows from Eq. 

d 2 -i d 2 ~i 
£\E(i)) = J2 S a( 3\E a )) «-> S(Ep) = Yl S *pE a , 

a=0 a=0 

it is easy to confirm that \ can be written as 

d 2 -l d 2 -l 

X= Y Xaf3\E Q ))((E p \= J2 S a0 E a ®E* . (20) 

a, 13=0 «:/3=0 



Equations l|19|) and l|20|l are one of the main results of this 

paper. From these equations, we find that S and x are 
complementary to each other in the sense that they can 
be interchanged if we exchange the two operator bases 
\E a )){{Efs\ and E a ® E* p on Hf 2 in their expressions. 
These equations show that the 5-matrix (^-matrix) is 

given by the expansion coefficients of S (x) with respect 
to \E a ))((Ep\ as well as those of x (S) with respect to 
E a © E%, which is explicitly written as 

Xa! 3 = Tr(|^ a ))«^|)t| = Tr(£ Q © E^S, (21) 

S aP = Tr(|E a ))«^|)tj = T r (E a © (22) 

From Eqs. I|19|) and i|20[l . we can explore the mutual 
conversion formulas between x~ and 5-matrices. To this 
end, let us define a bijection between the two operators 
on Hf 2 originally found by Havel [iof . 

d 2 -i 

A (°) = H (i®n y )@(ity®i), (23) 

7=0 

which is also considered to be the super-superoperator 
acting on TLd- Then, we have the following Theorem [40| . 

Theorem 2 (Havel (2003)) For arbitrary operators X 
and Y in TLSd, we have 

|1})«Y|=A(1©Y*), (24) 

X®Y* = h{\X)){{Y\). (25) 

Theorem |2 connects two relevant operators on TLf 2 hav- 
ing different tensor product structures, i.e., the Kro- 
necker product X © Y* and the dyadic product |X))((Y|. 
To prove the Theorem [21 we first note the following 
lemma. 

Lemma 3 The identity operator on TL® 2 and the (un- 
normalized) density operator of the isotropic state on 
TL® 2 are related as follows. 

|7))«J|=A(J®J), (26) 

/®/ = A(|/))«J|). (27) 

It is straightforward to confirm Eqs. (|26|l and l|27|l by 

writing and I © I using the standard state basis 

\i) explicitly. Then, it follows that 

|X»«Y| = (X®I)\I)){{I\(I®Y*) 

= A((X © © © Y*)) = A(X © Y*), 
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x®y* = (X(g>i)(i(g>i)(i(g>Y*) 

= A((X f)|7))«7|(/ f •)) = A(|l)(y|), 

where we used Eq. (j^J. Accordingly, Theorem is 
proved. At this point, we note that the action of the 
bijection A (0) corresponds to reshuffling of the matrix 
introduced by Zyczkowski and Bengtsson: if we con- 
sider the matrix for the operator on Ti.f 2 defined with 
respect to the standard state basis, the mapped opera- 
tor by A (0) has a reshuffled one of the original matrix 
pfij . The bijection A(Q) is also closely related to the ma- 
trix realignment introduced by Chen and Wu to discuss 
the separability criterion for the bipartite density matrix 
[3^| . It is easy to confirm that A(©) is involutory, that is, 
A(A(0)) = 0. It should be also noted that A(0) does not 
preserve Hermiticity and the rank of the transformed op- 
erator, so its spectrum is not preserved [2(j. This means 
that A(©) represents a non-physical operation. It follows 

from Theorem[2]that x — A(«S) and S = A(x), i.e., S and 
X are bijective. 

From this bijective relation, we can explore the bijec- 
tion between x~ and 5-matrices. To this end, we expand 
\E a ))((Ep\ in terms of E a ®E*p, and vice versa. Since 
they are bijective, it follows that 

E a ®E* = A(\E a }}{{E \) 
d 2 -i 

= Yl \Ea>))((Efi>\Ma>(y, a (i, (28) 

a',/3', 7= 

\E a ))((Ep\ = A(E a ®E*p) 

= E a , ® E* pl {M^) a 

'0',a,0, ( 29 ) 

a',/3'=0 

where M a ipi j0l p is a! ft; a/3-entry of the d 4 x d 4 complex 
matrix M. It is explicitly given as 

Tr(\E a ,}}{{E ,\^E a ®E*p 
d 2 -i 

= E^'A" (30) 

7=0 

where the coefficients 

Qlp = {{E a \{I®^)\Ep)), (31) 



Rl p = {{E a \{%-,®i)\E fj )) (32) 

are the computable matrix elements of the operators 

1 7r 7 and 7r 7 1 on Tif 2 defined with respect to 
,2 -. 

the basis set {\E a ))}'^S . These coefficients form 
the d 2 x d 2 complex matrices Q 1 = [Q^]^ g~= and 
FT = [i2^o]^ gl Q . It is easy to confirm that M is Hermi- 
tian as well as unitary, i.e., (M^) a /p/ t0l p = {M a p tCt /p')* — 



Ma'p'^p and J2 a " ,p»( M ^) a p,a"p>> M a »p„. a ,p, = 
J2 a ",p" M aP;a"p"(M^) all p,^ a ,p, = 5 aa >Spp,. By us- 
ing these results, we obtain a bijection between the x~ 
and S -matrices: 

d 2 -i 

X= £ Q 7 SR\ (33) 

7=0 

d 2 -i 

5=E Q J XR\ (34) 
7=0 

which are given by the sum of the multiplication of three 
known matrices Q 7 , x, and i? 7 , and evidently com- 
putable. 

The above formulation can be straightforwardly ex- 
tended to describe two-qudit operations. In this case, 
different choices of basis sets are allowed for two sys- 

d 2 — 1 

terns. Let us choose the set {E a } a=a acting on the first 

d 2 — 1 

qudit space and the set {F a } a=0 acting on the second 
qudit space. The general two-qudit superoperator can be 
written as 

S= J2 s «P,^ a ))\F0)){0 7 \{{Fs\ (35) 

a,/3,7,<5=0 

acting on the d -dimensional L-space L^i as well as 
d 2 -i 

5(0)= X a p^sE a ®FpQE\®Fl (36) 

a, /3 : 7, (5=0 

acting on HS d 2. These superoperators are character- 

d 2 — 1 

ized with the d x d 4 matrices S = [S a p^$] a „ S=Q and 

d 2 -l 

X = [Xa0,js] at p^ tS=O ' The bijective Choi operator is de- 
fined on the <i 4 -dimensional -ff-space Tif 4 that is identi- 
fied with the L-space as follows: 

.(13) s (24) 

X ^S 0X {d 2 pf 2) ®pf\ (37) 

where the indices refer to the factors in 7i® in which the 
corresponding operations have a nontrivial action |4lll42l 

IZ3I ] . Then, it is straightforward to show that S and x can 
be written as follows: 

d 2 ~\ 

S = £ S Q p n s\E a ))\Fp))((E,\((F s \ 

a,/3,7,5=0 
d 2 -l 

= X«0,isE a ®E*®Fp®FI. (38) 

a, P, 7, 8=0 
d 2 -l 

X = £ X «w\E«))\Fp)}{{E^{{F s \ 

a,/3,7,<5=0 
d 2 -l 

= £ S a p^ s E a ®E*®Fp®F£, (39) 

a,/3,7,<5=0 



7 



It follows that the two operators are certainly bijective: 

X = A ® A(<S) and S =A ® A(x). From these bijective 
relations, we can explore the bijection between the x~ 
and 5-matrices for a two-qudit operation as 

d 2 -i 

X= Q J ®S X SR~< ®T X , (40) 

7,A=0 
d 2 -l 

S= Q~< ®S X X R~< ®T X , (41) 

7,A=0 

where d 2 x d 2 matrices Q 1 and PC are given by Eqs. 
ijjHJ and Q32JI . and the matrix entries of S x and T A are 
given by the matrix elements of the operators I ® 7r 7 
and 7r 7 (8 I on Hf 2 defined with respect to the basis set 

SZp = {{F a \{i®^)\Fp)), (42) 

^ = ((4|(^®/)|^)). (43) 

We can further extend the above formulation to describe 
n-qudit operations. In this case, the bijective relation 

between S and x reads x = A® ra (5) and S =A®"(x). By 
using these relations, we can straightforwardly extend 
Eqs. H4Ufl an d <|41|) to the bijection between the x~ ancl 
5-matrices for an n-qudit operation. 

It is obvious that not all the space of x~ an d 5-matrices 
corresponds to physically realizable operations. For ex- 
ample, we can describe an anti-unitary operation by us- 
ing the x~ an d 5-matrices, which is evidently an unphysi- 
cal operation. The requirement for the x~ an d 5-matrices 
to represent physical quantum operations has been ex- 
tensively studied by many researchers. In the following, 
the requirements common for the single- and two-qudit 
operations are summarized |2l| . 

Condition 1 (Hermiticity) The physical quantum oper- 
ation 6 should preserve Hermiticity; i.e., & maps any 
Hermite operator into an Hermite operator. 

Condition 2 (Positivity) The physical quantum opera- 
tion 6 should be positive; i.e., & maps any positive op- 
erator into an positive operator. 

Condition 3 (Complete positivity) The physical quan- 
tum operation & should be completely positive; i.e., posi- 
tivity is preserved if we extend the L-space and HS-space 

by adding more qudits. That is, the superoperator S (g> I 
on the extended spaces should be positive. 

It is known that is sufficient for ^ and [21 and [21 is suffi- 
cient for^ Therefore, we require complete positivity for 
a physical quantum operation. Complete positivity can 
be expressed as a particularly simple condition for the 
X-matrix. 



Theorem 3 The linear operation & is completely posi- 
tive, iff the x-matrix is positive. 

This is natural on physical grounds, because the Choi 
operator x should be an unnormalized density operator 
associated with the system which was subjected to the 
quantum operation as will be discussed in the next sec- 
tion. In addition to [31 any physical quantum operation 
should satisfy the following condition. 

Condition 4 (Trace non-increasing) The physical quan- 
tum operation 6 should be trace non-increasing; i.e., the 
mapped operator should have trace less than one. 

This condition is simply expressed as the restriction on 
the ^-matrix: Tiix < I^ 1 ' for a single-qudit operation 
and Tri3X < 1^ f° r a two-qudit operation, where j' n J 
is an identity matrix with size d n . 

It is needless to say that the x~ an d 5-matrices can 
be defined with respect to arbitrary operator basis sets. 
Once these matrices are given with respect to a particu- 
lar operator basis set, they can be converted into those 
defined with respect to the other basis set. It is obvious 
from Eqs. flty. l(2UJl. and that the two matri- 
ces defined with respect to different bases are unitarily 
equivalent. To be specific, let x E an d S E be the x~ an d 
5-matrices for a single-qudit operation defined with re- 
spect to the operator basis set {E a }^~Q , and x F an d 
S F be those defined with respect to the operator basis 
set {F a } a=0 , where the two bases are unitarily related 
as shown in Eqs. H15fl an d Then, these matrices 

should be written as 

s f = u \ S e u ^ ( 44 ) 



X F =U^x E U, (45) 

where U = [l^ap]^^ is a eP x d 2 unitary matrix. For the 
case of two-qudit operation, we need to extend the set 
of the operator basis to cover all the possible basis sets 
defined for the two-qudit operator space, that is, one that 
is a factorable set as well as not a factrorable set with 
respect to the first and second systems. The general set 
of operator basis {<I > 7 } 7= ^) 1 on Hf 2 should be unitarily 
related to the factorable operator bases E a ® Fp. If we 

1 

introduce a d 4 x d 4 unitary matrix U = [W a /j] tt=0 that 
relates l> 7 and E a ® Fp: 

d 2 ~i 

*7= E a ®F p U [aAl , (46) 

a, 0=0 

where [a, (3] := d 2 a + [3, it follows that Eqs. and 
(14511 also hold for two-qudit operations. 

The x~ an d 5-matrices can be diagonalized by choosing 
the appropriate operator basis sets, but are not neces- 
sarily diagonalized simultaneously by a unique set. The 
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operator basis set that diagonalizes the x-matrix, each 
element of which is multiplied by the square root of the 
associated eigenvalue, forms a particular set of Kraus op- 
erators in the Kraus form of the quantum operation. Any 
set of Kraus operators can be obtained by noting the uni- 
tary freedom in the Kraus form [2|. It follows from Eqs. 
inH, 1201), jnH: and (J2HI that the same operator basis 
set with the associated set of eigenvalues also gives an 
operator-Schmidt decomposition for the L-space super- 
operator S [23|. Therefore, the Kraus rank for the HS- 
space superoperator S (©) and Schmidt number of the 
L-space superoperator S must be equal. 



IV. APPLICATIONS OF X - AND S-MATRICES 

This section presents the several applications of the \- 
and 5-matrices. We discuss how these matrices and the 
present formulation are useful for analysis and design of 
quantum operations. 

Experimental identification of quantum opera- 
tions 

In the first example, we explain how useful the present 
formulation is for experimental identification of quantum 
operations % 0, H |l E3, El 03 This task is impor- 
tant because the development of any quantum device 
or circuit for quantum computation and communication, 
which can be considered as an input-output system that 
performs an intended quantum operation on its input 
state and transforms it into its output state, necessarily 
requires experimental benchmarking of its performance. 
The identification of a two-qudit device is particularly 
interesting from a practical viewpoint as well as a scien- 
tific one because it may involve a nonseparable operation 
which has a purely quantum mechanical nature, i.e., it 
cannot be simulated by using any classical method. 

Identification of an input-output system amounts to 
identifying its \~ or 5-matrix, since these matrices char- 
acterize the system in question completely as far as input 
and output data are concerned. The evaluated matrices 
should reproduce the behavior of the system well enough 
when the system is stimulated by any class of inputs of 
interest, and they should be useful for engineering the 
system of interest, e.g., to permit control of the system, 
to allow transmission of information through the system, 
to yield predictions of future behavior, etc. Identification 
problems are commonly regarded as inversion problems, 
where the x~ 01 5-matrix is to be statistically estimated 
from incomplete prior knowledge of the system, using 
prior knowledge of corresponding inputs and the collec- 
tion of data obtained by measurement of outputs that 
usually contain noise. In what follows, it will be shown 
that this common belief is not the case for identification 
of quantum operations. To be specific, we can estimate 
both the x~ an d 5-matrices without any inversion pro- 
cedure if we can make use of an entangled resource and 



a sequence of local measurements assisted by classical 
communication 0H [H IH El • 

Consider first the identification of a single-qudit opera- 
tion. Equations i|I9|) and i|20|) show that all the elements 
of the x~ an d 5-matrices are given by the expansion co- 
efficients of the operators x or 5 with respect to two 
different types of operator basis on Tif 2 - Of these two 
operators, the Choi operator x is particularly useful since 
it is a positive operator associated with the physical state 
of the bipartite object. To be specific, it can be inter- 
preted as the unnormalizcd output state from the system 
in question where qudit 1 of the two qudits prepared in 
an isotropic state is input into the system and undergoes 
the quantum operation 6 while qudit 2 is left untouched. 
Therefore, we can prepare the output state corresponding 

to the normalized Choi operator x/d = S®I(pi) with the 
use of several copies of the isotropic-state input for the 
two qudits. Thus, the identification of a single-qudit op- 
eration reduces to the identification of a two-qudit state. 
It follows from Eq. 11' 21) that every element of the iS- 
matrix can be directly obtained by determining the ex- 
pectation value of the corresponding product operator 
basis (E a <g> Eg) for the output states after the quantum 

operation has taken place. If the basis E a is chosen to 
be the Hermitian operator basis A Q , it suffices to make 
a set of d 4 -independent local measurements assisted by 
classical communication to determine the whole set of the 
real expectation values (A a £g> A^). Accordingly, we can 
obtain the 5-matrix defined with respect to the Hermi- 

d 2 — l 

tian operator basis set {A Q } a=0 . Once the 5-matrix is 
obtained, it is easy to convert it to the x-matrix defined 
with respect to the same basis by using Eqs. H3III - 1I33|) 
and also into the x~ an d 5-matrices defined with respect 
to the arbitrary chosen basis by using Eqs. i|44[l and J45I . 

The identification of a two-qudit operation can be 
carried out in the same way as the identification of a 
single-qudit operation. In this case, we prepare the state 
corresponding to the Choi operator in Eq. (|37[l with the 
use of several copies of the product of isotropic states 
prepared in the four qudits. To prepare the output 
state, we initially prepare the product of isotropic states 
p\ 12 ^ pf 4 " 1 in two pairs of two qudits (qudits 1-2 and 
qudits 3-4). Then qudit I and qudit 3 are input into the 
system in question, undergo the quantum operation 6 

jointly while the other qudits arc left untouched. This 

- 2 . 

setup leads to an output state x/d in four-qudits. Thus, 
the identification of a two-qudit operation reduces to the 
identification of a four-qudit state. It follows from Eq. 
(|38|l that every element of the 5-matrix for the two-qudit 
operation can be directly obtained by determining the 
real expectation value of the corresponding product 
operator basis (A Q ® A 7 <g> A^ <g> X5) for the output states 
after the quantum operation has taken place, if all the 
relevant basis sets are chosen to be the Hermitian opera- 
tor basis A Q . It suffices to make a set of <i 8 -independent 
local measurements assisted by classical communication 
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to determine the whole set of real expectation values 
(A Q <g> A 7 g> X/3 <g> Xs). Accordingly, we can obtain the 
5-matrix defined with respect to the Hermitian operator 

basis {A a C§> \p} a « =0 - The S-matrix can be converted 
into the x- matr i x by using Eqs. (|4*Ull - (l4"3)l . The S- 
and %-niatrices defined with respect to arbitrary chosen 
bases can be obtained by applying the appropriate 
matrix unitary transformation. 

Matrix analysis of quantum operations 

This section discusses in what way the x~ an( i $~ 
matrices contribute to developing quantum devices and 
circuits for quantum computation and communication. 
We consider two classes of applications in which these 
matrices offer useful mathematical models for quantum 
operations. The first one concerns physical and infor- 
mation theoretic analysis of quantum operations and the 
other concerns a logical calculus of quantum circuits or 
algorithms comprised of a sequence of quantum opera- 
tions. 

Let us first consider the physical and information the- 
oretic analysis of quantum operations. For this purpose, 
it is preferable to use the ^-matrix. This stems partly 
from the fact that the ^-matrix is positive and isomor- 
phic to the density matrix in the doubled Hilbert space. 
Physically, the diagonal elements of the process matrix 
show the populations of, and its off-diagonal elements 
show the coherences between, the basis operators mak- 
ing up the quantum operation, analogous to the inter- 
pretation of density matrix elements as populations of, 
and coherences between, basis states. Owing to Jami- 
olkowski isomorphism, the dynamic problems concerning 
quantum operations can be turned into kinematic prob- 
lems concerning quantum states in a higher dimensional 
space, and one can make use of a well-understood state- 
based technique for analyzing the quantum operation. In 
what follows, we show several illustrative examples and 
interesting problems from the physical and information 
theoretic viewpoints. 

The first example concerns the fidelity or distance mea- 
sure between two quantum operations. Several measures 
that make use of the above isomorphism have been pro- 
posed to quantify how close the quantum operation in 
question is to the ideal operation (usually a unitary op- 
eration) we are trying to implement. For example, the 
state fidelity defined between two states is extended to 
compare the two operations. The process fidelity F p is 
defined by using the ^-matrix \ of the system in ques- 
tion and the rank one ^-matrix Xideai of the ideal system 
in the state-fidelity formula, that is, F p = -^TrxXideai, 
where n~l for single-qudit operation and n=2 for a two- 
qudit operation. The average gate fidelity F defined as 
the state fidelity between the output state after the quan- 
tum operation and the ideal output can be calculated 
from the process fidelity. The purity defined for the den- 
sity matrix can be extended to characterize how much 
of a mixture the quantum operation introduces, which is 
also represented by the simple function of the %-matrix 



of the system in question [HI 03 EH ■ 

The next example concerns the analysis of a quantum 
operation acting on the composite system. As mentioned 
before, the ^-matrix of a two-qudit operation is inter- 
esting from a practical as well as a scientific viewpoint. 
The Jamiolkowski isomorphism for a two-qudit opera- 
tion (Eqs. (|S7|) and implies that the notion of 
entanglement can be extended from quantum states to 
quantum operations. Analogously to what happens for 
states, quantum operations on a composite system can 
be entangled |46l |47| . A quantum operation acting on 
two subsystems is said to be separable if its action can 
be expressed in the Kraus form 

J(0) =^(i i ®A)0(ii®A) t , (47) 



where A, and Bi are operators acting on each subsystem 
|48l I49L IHcj. Otherwise, we say that it is nonseparable 
(or entangled). Quantum operations that can be per- 
formed by local operations and classical communications 
(the class of LOCC operations) are described by separa- 
ble quantum operations, yet there are separable quantum 
operations that cannot be implemented with LOCC op- 
erations with probability one |4lL l43l l51j . Anyway, these 
are useless for creating entanglement in an initially un- 
cntangled system. It has been pointed out by several 
authors that the separability and entangling properties 
of quantum operations acting on two systems can be dis- 
cussed in terms of the Cho i op erator for two-qudit op- 
erations (Eq. ((SHI) In the present context, 
this reduces to discussing the separability properties of 
the ^-matrices. For example, there is a condition for 
the x-matrix equivalent to Eq. 1471) : a quantum opera- 
tion acting on two subsystems is separable if its %-matrix 
can be written as x = TliixY^ 1 ® x\ B \ where Xi an d 
x\ are the ^-matrices for the quantum operation act- 
ing on each subsystem. Thus, the separability of general 
quantum operations acting on the composite system is 
reduced to the separability of its x-matrix. Since the 
separability criterion and measure for the general d A x d 4 
positive matrix is not fully understood, it remains as an 
important problem for quantum information science to 
find such a criterion and measure for general two-qudit 
quantum operations. 

Let us turn our attention to a logical calculus of quan- 
tum circuits or algorithms. For this purpose, the 5- 
matrix is practically useful. This follows from the fact 
that i^space superoperator algebra works just like Dirac 
operator algebra. For example, consider the scenario in 
which two-quantum operations ©1 and 62 act sequen- 
tially on a quantum system. Assume that the associated 
5-matrices are given with respect to the same operator 
basis set {\E a ))}a=Q in Cn, where N = d 2 for single- 
qudit operation and N = d for two-qudit operation. 
Then the composite operation & = &± o 6 2 is described 
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by the multiplication of Z^space superoperators 

AT-l 

i = JiJ 2 = ^ 5^|E Q ))((^|, (48) 

q,/3=0 

where 

JV-l 

S a p = Y, ( 5 i) Q7 (<5 2 ) 7/3 . (49) 

7=0 

The extension to the case in which a sequence of a finite 
number of quantum operations is applied to the same 
quantum system is trivial. Equation (|49|l implies that 
the iS-matrix of the composite operation reduces to the 
multiplication of the 5-matrices of the individual opera- 
tions. This makes it practically advantageous using the 
5-matrix to make a logical calculus of quantum circuits or 
algorithms comprised of a sequence of elementary single- 
and two-qudit quantum operations. 

Consider next the quantum circuit or algorithm com- 
prised of a sequence of quantum operations each of which 
acts not necessarily on the same quantum system. This 
offers a general model for the quantum circuit acting on 
the large numbers of qudits To analyze and design 
such a quantum circuit, we need to consider quantum 
operations acting on the whole set of qudits and associ- 
ated extended 5-matrix. Such an extended 5-matrix is 
non-trivial, but its bijective ^-matrix is trivially obtained 
by taking the tensor product with identity matrix, that 
is, the x-matrix for the identity operation on the irrele- 
vant system. If the %-matrix for the quantum operation 
& on the relevant space is given by %, the x-matrix for 
the extended quantum operation is given by \ ® !• O n 
the other hand, we can trivially extend the conversion 
formulas (|40|l and l|41|) to those formulas for the quan- 
tum operation acting on more qudits. Therefore, we can 
calculate the extended 5-matrix for each quantum op- 
eration from the associated 5-matrix for the quantum 
operation acting on the relevant space. The 5-matrix for 
a sequence of operations acting on the space of the whole 
quantum systems can be calculated by multiplication of 
the extended 5-matrices for the individual quantum op- 
erations. 

The 5-matrix analysis of the quantum operation has 
the following potential advantage. It can deal with a 
non-unitary operation in which mixed state evolution 
occurs. Noisy quantum operation, probabilistic subrou- 
tines, measurements, and even trace-decreasing quantum 
filters can be treated. This is in contrast to the usual 
analysis based on unitary matrix which can deal only 
with unitary gate in which only the pure state evolution 
is allowed. It thus offers us an mathematical model to 
analyze and design the logical operation of wider range 
of the complex quantum circuits and algorithms . 

In the above discussion, we considered two applica- 
tional classes, i.e., physical and information theoretic 
analysis and logical calculus of quantum operations, in 
which the x~ an d 5-matrices matrices offer useful mathe- 
matical models. They have their own useful applications. 



The present formulation will offer us the way of building 
bridges across the two applicational classes. For exam- 
ple, the entangling properties of the quantum circuits 
comprised of a sequence of single- and two-qudit quan- 
tum operations acting on several qudits can be discussed. 
The present formulation will also help us to analyze and 
benchmark the quantum operation realized in the actual 
device. 



V. CONCLUSIONS 

We have considered two matrix representations of 
single- and two-qudit quantum operations defined with 
respect to an arbitrary operator basis, i.e., the x~ 
and 5-matrices. We have provided various change-of- 
representation formulas for these matrices including bi- 
jections between the %- and 5-matrices. These matrices 
arc defined with the expansion coefficients of two oper- 
ators on a doubled Hilbert space, that is, the Z^space 
superoperator and the Choi operator. These operators 
are mutually convertible through a particular bijection 
by which the Kronecker products of the relevant opera- 
tor basis and the dyadic products of the associated state 
basis are mutually converted. From this fact, the mutual 
conversion formulas between two matrices are established 
as computable matrix multiplication formulas. Extention 
to multi-qudit quantum operation is also trivial. These 
matrices are useful for their own particular classes of ap- 
plications, which might be interesting from a practical as 
well as a scientific point of view. 

We have presented possible applications of the present 
formulation. By using the present formulation, an exper- 
imental identification of a quantum operation can be re- 
duced to determining the expectation values of a Hermi- 
tian operator basis set on a doubled Hilbert space. This 
can be done if we prepare several copies of the isotropic- 
state input or the product of isotropic states input. By 
using the x-matrix, we can make a physical as well as 
a quantum information theoretic characterization of the 
quantum operation. In particular, the x-matrix is useful 
to discuss the entangling properties of the quantum oper- 
ation acting on the composite system, since the problem 
of the separability of the quantum operation is reduced 
to the problem of the separability of the x-matrix. On 
the other hand, the iS-matrix is useful when we discuss 
the typical quantum circuit comprised of a sequence of 
single- and two-qudit quantum operations each of which 
acts on different quantum qubits. It is possible by consid- 
ering the extended 5-matrix of each quantum operation 
acting on the whole state space of the relevant qudits. 
Such extended 5-matrices can be calculated from the as- 
sociated, bijective x-matrices by taking the tensor prod- 
uct with the appropriate identity matrix. Accordingly, 
we can calculate the 5-matrix for a quantum circuit by 
multiplying the extended 5-matrices of each operation. 
This should be very useful to analyze and design a wide 
range of the quantum circuits and algorithms involving 



11 



non-unitary operation. the CREST program of the Japan Science and Tcchnol- 

We thank Satoshi Ishizaka and Akihisa Tomita for ogy Agency, 
their helpful discussions. This work was supported by 



[1] K. Kraus, States Effects and Operations: Fundamental 
Notions of Quantum Theory, Lecture Notes in Physics 
Vol. 190 (Springer- Ver lag, Berlin, 1983). 

[2] M.A. Nielsen, I.L. Chuang, Quantum Computation and 
Quantum Information (Cambridge University Press, Ox- 
ford, 2000). 

[3] D. Aharonov, A. Kitaev, and N. Nisan, e-print 

quant-ph/9806029 
[4] P. Stelmachovic and V. Buzek, Phys. Rev. A 64. 062106 

(2001); 67, 029902(E) (2003). 
[5] H. Hayashi, G. Kimura, and Y. Ota, Phys. Rev. A 67, 

062109 (2003). 

[6] Y. Nambu et al., Proc. SPIE Int. Soc. Opt. Eng. 4917, 
13 (2002). 

[7] J. B. Altepeter et al., Phys. Rev. Lett. 90, 193601 (2003). 
[8] M. W. Mitchell, C. W. Ellenor, S. Schneider, and A. M. 

Steinberg , Phys. Rev. Lett. 91, 120402 (2003). 
[9] F. De Martini et al., Phys. Rev. A 67, 062307 (2003). 
[10] J. L. O'Brien et al., Phys. Rev. Lett. 93, 080502 (2004). 
[11] V. Secondi, F. Sciarrino, and F. De Martini, Phys. Rev. 

A 70, 040301 (R) (2004). 
[12] Y. Nambu and K. Nakamura, Phys. Rev. Lett. 94, 010404 
(2005). 

[13] C. M. Caves, J. Superconductivity 12, 707 (1999). 

[14] V.E. Tarasov, J. Phys. A. 35, 5207 (2002). 

[15] L. P. Hughston, R. Jozsa, and W. K. Wootters, Phys. 

Lett. 183, 14 (1993). 
[16] G. M. D'Ariano, P. Lo Presti, and M. F. Sacchi , Phys. 

Lett. A 272, 32 (2000). 
[17] I. L. Chuang and M. A. Nielsen, J. Mod. Opt. 44, 2455 

(1997). 

[18] A. G. White et.al., e-print |qua"nt-ph /0308lT5[ 

[19] E. C. G. Sudarshan, P. M. Mathews, and J. Rau., Phys. 

Rev.,121, 920 (1961). 
[20] K. Zyczkowski and I. Bengtsson, Open Syst. Inf. Dyn. 

11, 3-42 (2004). 
[21] P. Arrighi and C. Patricot, Annals Phys. 311, 26 (2004). 
[22] K. Blum, Density Matrix Theory and Applications 

(Plenum Press, New York, 1981). 
[23] A. Royer, Phys. Rev. A 43, 44 (1991); 45, 793 (1992). 
[24] S. W. Hawking, Phys. Rev. D 14, 2460 (1976). 
[25] R. M. Wald, Phys. Rev. D 21, 2742 (1980). 
[26] A. Fujiwara and P. Algoet, Phys. Rev. A 59, 3290 (1999). 
[27] A. Ben-Reuven, Phys. Rev. 141, 34 (1966) . 



[28] A. Ben-Reuven, Adv. Chem. Phys. 33, 235 (1975). 
[29] A. Ben-Reuven, Phys. Rev. A 22, 2572 (1980); 22, 2585 
(1980). 

[30] Y. Prior and A. Ben-Reuven, Phys. Rev. A 33, 2362 
(1986). 

[31] G. M. D'Ariano, and P. Lo Presti, Phys. Rev. Lett. 86, 

4195 (2001); 91, 47902 (2003). 
[32] M. Horodecki and P. Horodecki, Phys. Rev. A 59, 4206 

(1999). 

[33] B. M. Terhal and K. G. H. Vollbrecht, Phys. Rev. Lett. 
85, 2625-2628 (2000). 

[34] G. Mahler and V. A. Weberrufi, Quantum Networks: Dy- 
namics of Open Nanostructures, 2nd ed. (Springer, New 
York, 1998). 

[35] C. Miquel, J. P. Paz, and M. Saraceno, Phys. Rev. A 65, 
062309 (2002). 

[36] M. L. Aolita, I. Garcfa-Mata, and M. Saraceno, Phys. 
Rev. A 70, 062301 (2004). 

[37] M.-D. Choi, Linear Algebr. Appl. 10, 285 (1975). 

[38] A. Jamiolkowski, Linear transformations which preserve 
trace and positive semidefiniteness of operators, Rep. 
Math. Phys. 3, 275 (1972)! 

[39] K. Chen, and L-A.Wu, Phys. Lett. 306, 14 (2002); Quan- 
tum Info. Comp., 3 193 (2003). 

[40] T. F. Havel, J. Math. Phys. 44, 534 (2003). 

[41] J. I. Cirac, W. Diir, B. Kraus, and M. Lewenstein, Phys. 
Rev. Lett. 86, 544 (2001). 

[42] W. Diir and J. I. Cirac, Phys. Rev. A 64, 12317 (2001). 

[43] A. W. Harrow and M. A. Nielsen, Phys. Rev. A 68, 12308 
(2003). 

[44] M. Raginsky, Phys. Lett. A 290, 11 (2001). 

[45] A. Gilchrist, N. K. Langford, and M. A. Nielsen, e-print 

quant-ph/0408063 
[46] P. Zanardi, Phys. Rev. A 63, 040304(R) (2001). 
[47] X. Wang and P. Zanardi, Phys. Rev. A 66, 044303 

(2002). 

[48] V. Vedral, M.B. Plenio, M.A. Rippin, and P.L. Knight, 

Phys. Rev. Lett. 78, 2275 (1997). 
[49] H. Barnum, M.A. Nielsen, and B.W. Schumacher, Phys. 

Rev. A 57, 4153 (1998). 
[50] E. M. Rains, Phys. Rev. A 60, 173 (1999). 
[51] C.H. Bennett, D.P. DiVincenzo, C.A. Fuchs, T. Mor, 

E. Rains, P.W. Shor, J.A. Smolin, and W.K. Wootters, 

Phys. Rev. A 59, 1070 (1999). 



